WMAP 3yr data with the CCA: anomalous emission and impact of component separation on the CMB power spectrum

The Correlated Component Analysis (CCA) allows us to estimate how the different diffuse emissions mix in CMB experiments, exploiting also complementary information from other surveys. It is especially useful to deal with possible additional components. An application of CCA to WMAP maps assuming that only the canonical Galactic emissions are present, highlights the widespread presence of a spectrally flat "synchrotron" component, largely uncorrelated with the synchrotron template, suggesting that an additional foreground is indeed required. We have tested various spectral shapes for such component, namely a power law as expected if it is flat synchrotron, and two spectral shapes that may fit the spinning dust emission: a parabola in the logS - log(frequency) plane, and a grey body. Quality tests applied to the reconstructed CMB maps clearly disfavour two of the models. The CMB power spectra, estimated from CMB maps reconstructed exploiting the three surviving foreground models, are generally consistent with the WMAP ones, although at least one of them gives a significantly higher quadrupole moment than found by the WMAP team. Taking foreground modeling uncertainties into account, we find that the mean quadrupole amplitude for the three "good" models is less than 1 sigma below the expectation from the standard LambdaCDM model. Also the other reported deviations from model predictions are found not to be statistically significant, except for the excess power at l~40. We confirm the evidence for a marked North-South asymmetry in the large scale (l < 20) CMB anisotropies. We also present a first, albeit preliminary, all-sky map of the "anomalous" component.Comment: 14 pages, 17 figures, submitted to MNRAS, references adde

CMB signal in WMAP 3yr data with FastICA

We present an application of the fast Independent Component Analysis (FastICA) to the WMAP 3yr data with the goal of extracting the CMB signal. We evaluate the confidence of our results by means of Monte Carlo simulations including CMB, foreground contaminations and instrumental noise specific of each WMAP frequency band. We perform a complete analysis involving all or a subset of the WMAP channels in order to select the optimal combination for CMB extraction, using the frequency scaling of the reconstructed component as a figure of merit. We found that the combination KQVW provides the best CMB frequency scaling, indicating that the low frequency foreground contamination in Q, V and W bands is better traced by the emission in the K band. The CMB angular power spectrum is recovered up to the degree scale, it is consistent within errors for all WMAP channel combination considered, and in close agreement with the WMAP 3yr results. We perform a statistical analysis of the recovered CMB pattern, and confirm the sky asymmetry reported in several previous works with independent techniques.Comment: 10 pages, 7 figures, submitted to MNRA

Maximum likelihood, parametric component separation and CMB B-mode detection in suborbital experiments

We investigate the performance of the parametric Maximum Likelihood component separation method in the context of the CMB B-mode signal detection and its characterization by small-scale CMB suborbital experiments. We consider high-resolution (FWHM=8') balloon-borne and ground-based observatories mapping low dust-contrast sky areas of 400 and 1000 square degrees, in three frequency channels, 150, 250, 410 GHz, and 90, 150, 220 GHz, with sensitivity of order 1 to 10 micro-K per beam-size pixel. These are chosen to be representative of some of the proposed, next-generation, bolometric experiments. We study the residual foreground contributions left in the recovered CMB maps in the pixel and harmonic domain and discuss their impact on a determination of the tensor-to-scalar ratio, r. In particular, we find that the residuals derived from the simulated data of the considered balloon-borne observatories are sufficiently low not to be relevant for the B-mode science. However, the ground-based observatories are in need of some external information to permit satisfactory cleaning. We find that if such information is indeed available in the latter case, both the ground-based and balloon-borne experiments can detect the values of r as low as ~0.04 at 95% confidence level. The contribution of the foreground residuals to these limits is found to be then subdominant and these are driven by the statistical uncertainty due to CMB, including E-to-B leakage, and noise. We emphasize that reaching such levels will require a sufficient control of the level of systematic effects present in the data.Comment: 18 pages, 12 figures, 6 table

Blind component separation for polarized observations of the CMB

We present in this paper the PolEMICA (Polarized Expectation-Maximization Independent Component Analysis) algorithm which is an extension to polarization of the SMICA (Spectral Matching Independent Component Analysis) temperature multi-detectors multi-components (MD-MC) component separation method (Delabrouille et al. 2003). This algorithm allows us to estimate blindly in harmonic space multiple physical components from multi-detectors polarized sky maps. Assuming a linear noisy mixture of components we are able to reconstruct jointly the anisotropies electromagnetic spectra of the components for each mode T, E and B, as well as the temperature and polarization spatial power spectra, TT, EE, BB, TE, TB and EB for each of the physical components and for the noise on each of the detectors. PolEMICA is specially developed to estimate the CMB temperature and polarization power spectra from sky observations including both CMB and foreground emissions. This has been tested intensively using as a first approach full sky simulations of the Planck satellite polarized channels for a 14-months nominal mission assuming a simplified linear sky model including CMB, and optionally Galactic synchrotron emission and a Gaussian dust emission. Finally, we have applied our algorithm to more realistic Planck full sky simulations, including synchrotron, realistic dust and free-free emissions.Comment: 20 pages, 21 figures, 1 table, TeX file, accepted for publication in MNRA

Constraining Primordial Non-Gaussianity with High-Redshift Probes

We present an analysis of the constraints on the amplitude of primordial non-Gaussianity of local type described by the dimensionless parameter $f_{\rm NL}$. These constraints are set by the auto-correlation functions (ACFs) of two large scale structure probes, the radio sources from NRAO VLA Sky Survey (NVSS) and the quasar catalogue of Sloan Digital Sky Survey Release Six (SDSS DR6 QSOs), as well as by their cross-correlation functions (CCFs) with the cosmic microwave background (CMB) temperature map (Integrated Sachs-Wolfe effect). Several systematic effects that may affect the observational estimates of the ACFs and of the CCFs are investigated and conservatively accounted for. Our approach exploits the large-scale scale-dependence of the non-Gaussian halo bias. The derived constraints on {$f_{\rm NL}$} coming from the NVSS CCF and from the QSO ACF and CCF are weaker than those previously obtained from the NVSS ACF, but still consistent with them. Finally, we obtain the constraints on $f_{\rm NL}=53\pm25$ ($1\,\sigma$) and $f_{\rm NL}=58\pm24$ ($1\,\sigma$) from NVSS data and SDSS DR6 QSO data, respectively.Comment: 16 pages, 8 figures, 1 table, Accepted for publication on JCA

Estimating the tensor-to-scalar ratio and the effect of residual foreground contamination

We consider future balloon-borne and ground-based suborbital experiments designed to search for inflationary gravitational waves, and investigate the impact of residual foregrounds that remain in the estimated cosmic microwave background maps. This is achieved by propagating foreground modelling uncertainties from the component separation, under the assumption of a spatially uniform foreground frequency scaling, through to the power spectrum estimates, and up to measurement of the tensor to scalar ratio in the parameter estimation step. We characterize the error covariance due to subtracted foregrounds, and find it to be subdominant compared to instrumental noise and sample variance in our simulated data analysis. We model the unsubtracted residual foreground contribution using a two-parameter power law and show that marginalization over these foreground parameters is effective in accounting for a bias due to excess foreground power at low $\ell$. We conclude that, at least in the suborbital experimental setups we have simulated, foreground errors may be modeled and propagated up to parameter estimation with only a slight degradation of the target sensitivity of these experiments derived neglecting the presence of the foregrounds.Comment: 19 pages, 12 figures, accepted for publication in JCA

Component separation methods for the Planck mission

The Planck satellite will map the full sky at nine frequencies from 30 to 857 GHz. The CMB intensity and polarization that are its prime targets are contaminated by foreground emission. The goal of this paper is to compare proposed methods for separating CMB from foregrounds based on their different spectral and spatial characteristics, and to separate the foregrounds into components of different physical origin. A component separation challenge has been organized, based on a set of realistically complex simulations of sky emission. Several methods including those based on internal template subtraction, maximum entropy method, parametric method, spatial and harmonic cross correlation methods, and independent component analysis have been tested. Different methods proved to be effective in cleaning the CMB maps from foreground contamination, in reconstructing maps of diffuse Galactic emissions, and in detecting point sources and thermal Sunyaev-Zeldovich signals. The power spectrum of the residuals is, on the largest scales, four orders of magnitude lower than that of the input Galaxy power spectrum at the foreground minimum. The CMB power spectrum was accurately recovered up to the sixth acoustic peak. The point source detection limit reaches 100 mJy, and about 2300 clusters are detected via the thermal SZ effect on two thirds of the sky. We have found that no single method performs best for all scientific objectives. We foresee that the final component separation pipeline for Planck will involve a combination of methods and iterations between processing steps targeted at different objectives such as diffuse component separation, spectral estimation and compact source extraction.Comment: Matches version accepted by A&A. A version with high resolution figures is available at http://people.sissa.it/~leach/compsepcomp.pd

Planck Intermediate Results. IV. The XMM-Newton validation programme for new Planck galaxy clusters

We present the final results from the XMM-Newton validation follow-up of new Planck galaxy cluster candidates. We observed 15 new candidates, detected with signal-to-noise ratios between 4.0 and 6.1 in the 15.5-month nominal Planck survey. The candidates were selected using ancillary data flags derived from the ROSAT All Sky Survey (RASS) and Digitized Sky Survey all-sky maps, with the aim of pushing into the low SZ flux, high-z regime and testing RASS flags as indicators of candidate reliability. 14 new clusters were detected by XMM, including 2 double systems. Redshifts lie in the range 0.2 to 0.9, with 6 clusters at z>0.5. Estimated M500 range from 2.5 10^14 to 8 10^14 Msun. We discuss our results in the context of the full XMM validation programme, in which 51 new clusters have been detected. This includes 4 double and 2 triple systems, some of which are chance projections on the sky of clusters at different z. We find that association with a RASS-BSC source is a robust indicator of the reliability of a candidate, whereas association with a FSC source does not guarantee that the SZ candidate is a bona fide cluster. Nevertheless, most Planck clusters appear in RASS maps, with a significance greater than 2 sigma being a good indication that the candidate is a real cluster. The full sample gives a Planck sensitivity threshold of Y500 ~ 4 10^-4 arcmin^2, with indication for Malmquist bias in the YX-Y500 relation below this level. The corresponding mass threshold depends on z. Systems with M500 > 5 10^14 Msun at z > 0.5 are easily detectable with Planck. The newly-detected clusters follow the YX-Y500 relation derived from X-ray selected samples. Compared to X-ray selected clusters, the new SZ clusters have a lower X-ray luminosity on average for their mass. There is no indication of departure from standard self-similar evolution in the X-ray versus SZ scaling properties. (abridged)Comment: accepted by A&

Planck Intermediate Results. IX. Detection of the Galactic haze with Planck

Using precise full-sky observations from Planck, and applying several methods of component separation, we identify and characterize the emission from the Galactic "haze" at microwave wavelengths. The haze is a distinct component of diffuse Galactic emission, roughly centered on the Galactic centre, and extends to |b| ~35 deg in Galactic latitude and |l| ~15 deg in longitude. By combining the Planck data with observations from the WMAP we are able to determine the spectrum of this emission to high accuracy, unhindered by the large systematic biases present in previous analyses. The derived spectrum is consistent with power-law emission with a spectral index of -2.55 +/- 0.05, thus excluding free-free emission as the source and instead favouring hard-spectrum synchrotron radiation from an electron population with a spectrum (number density per energy) dN/dE ~ E^-2.1. At Galactic latitudes |b|<30 deg, the microwave haze morphology is consistent with that of the Fermi gamma-ray "haze" or "bubbles," indicating that we have a multi-wavelength view of a distinct component of our Galaxy. Given both the very hard spectrum and the extended nature of the emission, it is highly unlikely that the haze electrons result from supernova shocks in the Galactic disk. Instead, a new mechanism for cosmic-ray acceleration in the centre of our Galaxy is implied.Comment: 15 pages, 9 figures, submitted to Astronomy and Astrophysic

Planck intermediate results. III. The relation between galaxy cluster mass and Sunyaev-Zeldovich signal

We examine the relation between the galaxy cluster mass M and Sunyaev-Zeldovich (SZ) effect signal D_A^2 Y for a sample of 19 objects for which weak lensing (WL) mass measurements obtained from Subaru Telescope data are available in the literature. Hydrostatic X-ray masses are derived from XMM-Newton archive data and the SZ effect signal is measured from Planck all-sky survey data. We find an M_WL-D_A^2 Y relation that is consistent in slope and normalisation with previous determinations using weak lensing masses; however, there is a normalisation offset with respect to previous measures based on hydrostatic X-ray mass-proxy relations. We verify that our SZ effect measurements are in excellent agreement with previous determinations from Planck data. For the present sample, the hydrostatic X-ray masses at R_500 are on average ~ 20 per cent larger than the corresponding weak lensing masses, at odds with expectations. We show that the mass discrepancy is driven by a difference in mass concentration as measured by the two methods, and, for the present sample, the mass discrepancy and difference in mass concentration is especially large for disturbed systems. The mass discrepancy is also linked to the offset in centres used by the X-ray and weak lensing analyses, which again is most important in disturbed systems. We outline several approaches that are needed to help achieve convergence in cluster mass measurement with X-ray and weak lensing observations.Comment: 19 pages, 9 figures, matches accepted versio